Colloid aggregation induced by oppositely charged polyions 
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The "polymer reference interaction site model" (PRISM) integral equation formalism is used to 
determine the pair structure of binary colloidal dispersions involving large and small polyions of 
opposite charge. Two examples of such bidisperse suspensions are considered in detail, namely 
mixtures of charged spherical colloids and oppositely charged polyelectrolyte chains, and binary 
mixtures of oppositely charged large and small clay platelets. In both cases clear evidence is found 
for aggregation of the larger particles induced by the polyionic counterions, signalled by a strong 
enhancement of long wavelength concentration fluctuations. 



I. INTRODUCTION 

In metallurgy it has been known for several millenia 
that mechanical properties of many metals can be con- 
siderably improved by alloying two or several elements. 
Similar improvements of material properties of " soft mat- 
ter" can be achieved by mixing two or more macromolec- 
ular components in fluid solutions or dispersions. One 
of the oldest known examples is the stabilization of col- 
loidal suspensions by the addition of adsorbing polymer, 
which was already used in ancient Egypt for the prepara- 
tion of inks. On the contrary, flocculation of colloids can 
be driven by the additicm of non-adsorbing polymer via 
the depletion mechanismtl. More recent is the realization 
of the key role of Coulombic interactions in determin- 
ing and improving material or physiological properties 
of complex colloidal and biomolecular systems, and their 
exploitation in numerous technological and biomedical 
applications. Widely studied systems include suspen- 
sions of charged colloids and polyelectrolyteaa'cl, polymer- 
clay nanoepmpositesQ or DNA bundles and DNA-lipid 
complexesEl. For example, even small polyelectrolyte ad- 
ditives can have a substantial impact on the aggrejgation 
and kinetic stability of charged colloidal particleaJ. On 
the other hand, just a small fraction of claw, fillers doubles 
the tensile strength of a polymer matrixQ, while highly 
asymmetric mixtures of charged clay platelets are used 
to improve rheological properties of drilling fluidsQ. 

To design complex mixtures with optimal properties, 
it is essential to gain a quantitative understanding of 
correlations and effective interactions between various 
macromolecular species. This can be achieved using the 
well-established techniques of liquid-state theoryO. The 
PRISM (" polymer reference interaction site model" ) the- 
ory, originally designed for the study of polymer solutions 
and meltaJ, has been recently extended to investiftate 
various charged and uncharged colloidal systemsEHlEil, 
as well as, mixtures of neutral spherical colloids and 
polymeral3"E3 pid like-charged spherical colloids and 
polyelectrolytests. 

This paper focuses on the structure and stability of 
dispersions of charged colloidal particles in the presence 



of oppositely charged polyions, which play the role of 
mesoscopic counterions. There is a considerable litera- 
ture on the widely studied case of microscopic counteri- 
ons and the resulting effective interactions between the 
highly charged colloidal particlestZl. In particular, it is 
now well established that divalent counterions can lead 
to overscreening of the Coulombic interaction between 
spherical or rod- like colloidal particles, which results in a 
correlation-induced, short-range attraction between the 
latterllaO. The case where the counterions are meso- 
scopic or macromolecular polyions has received much less 
attention so far, except for the simplest case of polyelec- 
tolyte chains confined between two infinite parallel plates 
carrying charges of opposite sign to the polyelectrolyte, 
where overcompensation of the charge on the platelets, 
leading to effoetpe attraction between the latter, has 
been predictedEilO. 

Two systems involving mesoscopic counterions are ex- 
amined in the present paper, namely a) mixtures of 
highly charged spherical colloids and oppositely charged 
polyelectrolyte chains and b) binary mixtures of oppo- 
sitely charged platelets of widely different diameters. The 
key problem to be addressed is that of structural diag- 
nostics of a possible aggregation of the larger particles 
induced by the smaller particles, over a range of physi- 
cal conditions. This goal is achieved by determining pair 
correlation functions and partial structure factors within 
the framework of PRISM, which is adapted in this paper 
to the systems in hand. 



II. INTERACTION SITE MODELS 

The two systems under investigation are aqueous dis- 
persions or solutions, but in view of the mesoscopic scale 
of the particles, the solvent will be modelled as a struc- 
tureless dielectric continuum providing a macroscopic 
permittivity e. Any microscopic counterions or ions from 
added electrolyte will be considered at the linear response 
(or Debye-Hiickel) level, i.e., they will screen the electro- 
static potential due to the interaction sites on the colloids 
or polyions on a scale given by the usual Debye screen- 
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length AdE3. The underlying Born-Oppenheimer-like 
assumption entails that the charge distribution on the 
mesoscopic particles does not contribute to screeningc^. 
Since sites with charge of opposite sign are involved, the 
site-site interaction must contain a short-range repulsion, 
which will be modelled by hard sphere (HS) cores of di- 
ameter di. The interaction potential between sites on 
particles of species i and j, carrying the charges z^e and 
ZjC, will hence be of the generic form: 

f oo , r < Udi + dj) , 

^vir)^\ z,z,e^ . . u, , (2.1) 

t sr exp(-Ki3r) , r > ^(d, + dj) , 

where no = A^^. On the polyelectrolyte, Np charged 
sites are distributed at regular intervals along the lin- 
ear macromolecular backbone, each site being associated 
with a monomer or segment of diameter dp. The clay 
platelets will be modelled as circular discs of radius i?, 
carrying interaction sites uniformly distributed on the 
surface; the core diameter associated with these sites ac- 
counts for the finite thickness of the platelets. As regards 
the spherical colloids, two models may be considered. In 
model A, Ng charged sites are distributed over the spher- 
ical surface of radius R. For a sufhciently large density 
Ns/{^ttR^) of the surface sites, the hard cores of diameter 
ds associated with each site will ensure that the spheres 
are impenetrable. In model B, the impenetrability of 
the spherical colloids is ensured by an uncharged central 
(c) site, with a core radius identical to the particle ra- 
dius, while the electrostatic interactions are carried by 
Ns surface sites distributed on a spherical shell of radius 
R — ds/2. Each surface site carries a charge equal to 
the total charge on the colloid divided by the number of 
sites. Invoking Gauss' theorem, it might be argued that 
this distribution should be equivalent to a single site car- 
rying the full colloidal charge at the centre of the sphere. 
This equivalence would hold within an exact theory, but 
since an approximate closure relation is to be used, as 
discussed in section III, more accurate results may be 
expected when interaction sites with lower charges are 
used, such that the Coulomb coupling is weaker. 

Within the PRISM formalism to be discussed in the 
next section, all interaction sites of the same type (i.e., 
p, c, or s) are assumed to be equivalent, leading ±o a con- 
siderable simplification over the original RISME3. This 
assumption is obviously exact for spherical particles with 
surface sites only (model A), but is an approximation for 
polyelectrolytes of finite length (since end effects are ne- 
glected), and even more so for platelets, where edge ef- 
fects are not taken into consideration. However previous 
experience with neutral or charged rodalJ and plateletstil 
shows that the basic PRISM assumption does lead to sen- 
sible results for the pair structure of these rigid particles. 



III. MULTICOMPONENT PRISM 

Consider a multicomponent system involving v species 
of particles with number densities Pa ^ ot < v)^ each 
particle of a species involving Ua classes of exactly or ap- 
proximately equivalent interaction sites. The total num- 
ber of distinct classes of sites in the mixture is 



a=l 



= N. 



(3.1) 



It proves convenient to assign an index i (1 < i < N) to 
each class of sites, and to order them such that classes 
1 < i < rii belong to particles of species 1, classes rii+l < 
i < ni+n2 belong to species 2, ... , and classes N — n^ < 
i < N belong to species v. Each class i contains iV^ 
equivalent interaction sites, and if the positions of sites 
belonging to classes i and j are designated by = r[,^^ , 

and r; = rp^ < k < Ni\ I < I < Nj), then one may 
define the usual intramolecular structure factors, or form 
factors 
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(3.2) 



provided i and j sites belong to the same particle, while 
u!ij{q) = otherwise. The Wij(q) are hence elements of a 
N X N box-diagonal matrix uj(q). In the present paper 
rigid platelets and spherical colloids are considered, so 
that 
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where r^; = jr^ — r(|, g = |q|, and the prime signifies that 
the "self" term fc = / is to be left out for the diagonal 
elements i = j- PRISM is based on the assumption that 
all direct correlation functions between sites on pairs of 
different particles are identical if these sites belong to the 
same class on each of the two particles. This leaves a total 
of N{N -I- l)/2 independent direct correlation functions 
Cij (q) . Corresponding total correlation functions hij (q) 
are defined by averaging over the NiNj correlation func- 
tions between all pairs of equivalent sites belonging to 
classes i and j on two particles of the same or differ- 
ent species. The N x N matrices of Cij{q) or /iy (g) are 
symmetric, but it proves convenient to introduce auxil- 
iary non-symmetric matrices H(g'), C(g) and W(g), with 
elements 

Hij{q) = piNiNjhij{q) , Cij[q) = PiCij{q) , (3.4) 
W,j{q)^ ^/NJV]LU,jiq), (3.5) 

where pi is the number density of that species of parti- 
cles to which sites of category i belong. With these con- 
ventions, the set of N{N + l)/2 Ornsteiii-Z«i-mke (OZ) 
relations can be cast in the compact formt2lEaEa 
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H(g) = W{q)C{q) [W{q) + U{q)] (3.6) 

which may be solved for H{q) according to: 

H(q) = ([I - W(g)C(g)]-i - I) W{q) . (3.7) 

The set of N{N + l)/2 independent OZ relations must 
be supplemented by as many closure relations between 
each pair of total and direct correlation functions. For 



intersite distances r > dij = {di + dj)/2, we 
partially linearized Laria-Wu-Chandler closure 



?t the 
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1 + Xij (r) , 



X^J{r) > 0, 



(3.8) 



with 



Xij{r) = hnj{r) - LUii * (/3uy + Cy) * Wjj(r) , (3.9) 

where (3 = l/{kBT) is the inverse temperature. The su- 
perscript refers to the underlying reference model with 
hard core interactions only (i.e. , without the screened 
Coulomb interaction in Eq. (2.1)), while the asterisk * 
denotes a convolution product. The reference site-site 
correlation functions ^(r) and h'^Ar) are obtained by 
solving the set of OZ equations ( |3.6[ ) with the atomic 
Percus-Yevick closure for hard spheres, i.e. h'^j{r) = —1, 
r < dij] 4(r) = 0, r > dijU. The closure (3.8) takes 
the long-range asymptotics of the direct correlation func- 
tions Cij(r), in terms ofr-the molecular hypernetted chain 



(HNC) approximatiorocj, into account, but prevents the 
exponential increase of the total correlation functions 
hij{r) for distances correspondijis; to strong attraction 
between oppositely charged sitesc3. By construction, the 
total correlation functions and their first derivatives are 

continuous when Xii(^) = 0- 

The closed set of equations (^) - ( ^.9[ ) are solved nu- 
merically by a standard iterative procedure to obtain pair 



distribution functions gij{r) 
structure factors: 



1 -|- hij(r), and partial 



Sij{q) = ujij{q)6ij + pihij{q) . 



(3.10) 



IV. COLLOID-POLYELECTROLYTE MIXTURES 

In this binary mixture {v = 2), species 1, referred to 
with the index p, are semiflexible polyelectrolyte chains 
with Np equivalent segments carrying each a charge ZpC = 
Zpc/Np (where Z^e is the total charge of|-each chain). 
The form factor of the semiflexible chain iscJ: 



Op{q) = 1 + ]^ E 



X exp 



P n=l 

3 



(4.1) 



where dp is the segment length and Ip is the persis- 
tence length. The spherical colloids of radius Rc (species 
2), carry two categories of sites: Ns rigidly arranged 
surface sites of diameter ds carrying an electric charge 
ZgC — Zcc/Ns (where ZcC is the total charge of each col- 
loidal particle), and one central site, labelled with the 
index c, which is non-interacting (model A), or the cen- 
tre of a hard sphere repulsion of radius Rc (model B). In 
this paper only calculations based on model A will be pre- 
sented. The corresponding form factors are uJcc{q) = 1, 



(4.2) 
(4.3) 



where the double sums in Eq. ( |3.3[ ) have been replaced 
by double integrals corresponding to a continuous distri- 
bution of interaction sites on the surface of the sphere of 
radius Rc- Since this correspond s fo rmally to the limit 
A^s oo, the two form factors ( |4.2[ ) and (13) are ap- 




proximately related by LUss{q) 



^q). In practical cal- 



culations Ns = 60 was chosen. If no microscopic co- 
and counterions are present, charge neutrality relates the 
polyelectrolyte and colloid density via: 



ZcPc + ZpPp = . 



(4.4) 



In the absence of mic ros copic ions, = 0, and the po- 
tentials Uij (r) in Eq. ( ^.l| ) reduce to bare Coulomb poten- 
tials. With three classes of sites (p, s and c), the matrix 
W((7) is given by: 



LOpp{q) 

w((7) = ( iOss{q) ^cs{q) 
i^cs{q) i^cc{q) 



(4.5) 



There are a priori 6 independent direct and total corre- 
lation functions Cij{q) and hij{q) = p, s or c). How- 
ever, independently of any specific closure relation, the 
OZ equations lead directly to the following three relations 
between the total correlation functions: 



,(g) = uJcs{q)hcc{q) , 



hps{q) = UJcs{q)hpc{q) ■ 



(4.6) 
(4.7) 
(4.8) 



These relatiorLS reflect the fact that the "free-rotation" 
approximatioio for rigid molecules is in fact exact when 
the interaction sites are distributed on the surface of a 
sphere. Other decoupling schemes whiiekr-iiave recently 
been put forward for polymeric systemfH'Ea are also ex- 
act in that limit. 

We have considered mixtures of colloids of radius Rc — 
150 nm, and of oppositely charged semiflexible polyelec- 
trolyte chains of contour length L — 600 nm and per- 
sistence length Ip — 25 and 50 nm in water at room 
temperature [e^ / [ekBT) — 0.714 nm). 
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FIG. 1. Polyelectrolyte pair correlation function Ppp(r) 
for a mixture of semiflexible polyelectrolytes (contour length 
L = 600 nm) and oppositely charged spherical colloids (radius 
iic = 150 nm) for various valences Z-p = -Zc = 200, 400, 600 
(from bottom to top) at pc — 0.5 /im""^. The solid and dashed 
lines are calculated using the persistence length Ip = 100 nm 
and Ip — 50 nm, respectively. The inset displays the poly- 
electrolyte-colloid surface site correlation function gps{r) for 
Zc = -Z^ = 400 and L = 100 nm. 




150 



FIG. 2. Polyelectrolyte pair correlation function gpp{r) 
for a mixture of semiflexible polyelectrolytes (contour length 
L = 60 nm, persistence length Ip = 10 nm) and oppositely 
charged spherical colloids (radius Rc = 20 nm) for various 
valences Zp = —Zc = 50, 75, 100 (from bottom to top) at 
Pc = 5 X 10~^ nm~^. The inset displays gpp{r) for a mixture 
of charged dimers {L — 1 nm) and oppositely charged small 
spherical particles {Rc = 2 nm) fovr-Zc — —60 together with 
computer simulation data (squares)L3. 

The calculated pair distribution functions gpp{r) are 
shown in Fig. ^for three values of Zc- The interchain cor- 



relations are seen to increase dramatically with and 
to be much more pronounced for the more flexible chains. 
Qualitatively similar results hold for a mixture involving 
smaller colloidal particles and shorter chains, as shown in 
Fig. H The rehability of PRISM is illustrated in the inset 
of that Figure, where the PRISM predictiops for gpp{r) is 
compared to available Monte Carlo dataE3 for a mixture 
of charged spheres and oppositely charged dimers. The 
unusually large amplitude of the first peak in gpp{r) is the 
signature of an enhanced polyelectrolyte density near the 
oppositely charged spherical colloids, reminiscent n f t h e 
well-known " counterion condensation" of small iopscilEa. 
This condensation may drive phase separationEj which 
can be diagnosed by large fluctuations in the colloid den- 
sity. 

Examples of the colloid centre of mass structure fac- 
tor Sccio) are shown in Fig. ^ A central peak in the 
q — > limit is seen to build up as the colloid density pc 
increases. The enhancement of the density fluctuations, 
as characterized by Scdq = 0), is illustrated in Fig. |^, 
as a function of pc for several values of Z^. Although 
no divergence occurs, the tendency towards segregation 
is clear. In a naive random phase approximation (RPA) 
picture, values of Scc{Q) > 1 point to an eff'ective long 
range attraction between the colloidal particles. The ex- 
istence of an attraction between the colloidal particles, 
induced by the oppositely charged polyelectrolyte, can 
be more firmly established by a direct inversion of the 
pair distribution function gssir) of surface sites on the 
colloids. 
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FIG. 3. Colloid centre of mass structure factor Scc{q) 
for a mixture of semiflexible polyelectrolytes (L = 600 
nm, Ip — 50 nm) and oppositely charged spherical colloids 
{Rc — 150 nm, Zc — 1000, Zp = —800) for various densities 
Pc = 0.004, 0.04, 0.1, 0.2, 0.4 pm~''^ (from bottom to top). 
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FIG. 4. Colloid centre of mass structure factor at zero scat- 
tering vector Scc(O) as a function of colloid density, for a 
mixture of semiflexible polyelectrolytes (L = 600 nm) and 
oppositely charged spherical colloids {Rc = 150 nm) for vari- 
ous valences Zc — —Zp — 100, 200, 300 (from bottom to top). 
The solid and dashed lines are calculated using the persistence 
length Ip — 25 nm and Ip — 50 nm, respectively. 

In the limit pc ^ 0, this effective potential between 
surface sites would coincide with the potential of mean 
force, Vssir) = —ksT log gss{r), but at finite concentra- 
tion of colloids, an inversion based on an accurate closure 
relation must be used. In other words the inverse problem 
is solved, whereby the pair distribution function g^sir), 
as calculated from the solution of the multicomponent 
PRISM equations, is known, and the effective potential 
Vssir) is extracted from the one-component version of 
the theory, using the corresponding one-component OZ 
and closure relations. Note that the inversion does not 
require an iterative procedure. Similar inversions have re- 
cently been successfully used to determine effective inter- 
actions between the centres of mass of interacting poly- 
mer coilf£3, and— between the interaction sites of rodlike 
polyelect roly t esEZl. 

Examples of the effective potentials are shown Fig. ||. 
As the charge ZcC increases, an attractive well develops 
and deepens. Instability towards colloid clustering and 
genuine phase separation is expected to occur when the 
well-depth becomes of the order ksT, so that yet stronger 
Coulomb interaction, i.e., more highly charged polyelec- 
trolyte, would be needed to drive the phase separation. 
This is presently under investig atio n. A comparison with 
the bare Coulomb interaction (|]^) clearly demonstrates 
the enhanced screening effect of the polyelectrolyte as the 
charge increases. The same figure also shows the effective 
potential Vcdr) between the centres of the colloidal parti- 
cles, obtained from an HNC inversion of the centre to cen- 
tre pair distribution function gcc{f) for one value of Zc- 



-0.5 



11 ' 

1 1 

\ V ^ ^ ^ 
\ ^ ^ ^ ^ 


' 1 ' 1 ' 1 

- 


^/^^^'^^ 1 


2R =300 nm 

C 

I.I.I 



100 200 300 400 500 600 



r [nm] 

FIG. 5. Effective potential between surface sites on neigh- 
bouring colloids (solid line), and the bare Coulomb site-site 
interaction (dashed line) for a mixture of charged spherical 
colloids (radius Rc = 150 nm) in the presence of oppositely 
charged polyelectrolytes (contour length L — 600 nm, per- 
sistence length Ip = 50 nm) at pc — 5pm~'^. From top to 
bottom the solid lines correspond to increasing values of va- 
lences (Zc = -Zp = 200, 400, 600), while the dashed lines 
correspond to decreasing values of Zc- The dotted line repre- 
sents the effective potential Vcc{r) between the centres of the 
colloidal particles as derived from an HNC inversion of gcc{r) 
for Zc = -Zp = 200. 



V. BIDISPERSE CLAY SUSPENSIONS 

We now turn to binary mixtures of oppositely charged 
disc-shaped platelets of different radii, i?i and i?2 sus- 
pended in water, in the presence of microscopic co- and 
counterions. A'^i and N2 interaction sites are assumed to 
be uniformly distributed over the surfaces of the discs of 
species 1 and 2, and interact via the screened Coulomb 
potential ( |2.lD ; the hard core diameters di and d2 de- 
termine the effective thickness of the platelets. In the 
limit of a continuous distribution of sites, the form fac- 
tof&. iLji{q) = uni{q) (1 < i < 2) are easily calculated to 

where Ji denotes the cylindrical Bessel function of first 
order and all interaction sites are assumed to be equiv- 
alent, i.e., edge effects are neglected. This may be a 
more severe approximation for platelets than for spheres 
(where the equivalence is exact), or rigid or semiflexible 
polyelectrolyte chains, but our previous results for the 
one-component case show that it leads to good agree- 
ment with available simulation datatil. The correlation 
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function matrices tl{q) and C{q) are now 2x2 matrices, 
and W(q) is diagonal. 

We have solved the PRISM equations for platelet pa- 
rameters and physical conditions appropriate for a water 
base drilling mud, made up of large, negatively charged 
bentonite platelets, (index 6), and much smaller, pos- 
itively charged mixed metal hydroxide platelets (index 
m). The radii of the two species were taken to be 
Ri ~ Rb ^ lO'^ nm and i?2 = Rm = 25 nm; the ef- 
fective thicknesses were chosen to be db = dm = 10 nm, 
mimicking conditions where the platelets are not fully 
dispersed and form stacks. The strong screening in water 
base drilling muds is mainly due to the addition of NaOH 
(pH=12), and this is modelled by adopting an inverse 
screening length, kd = 1 nm"'^. A homogeneous interac- 
tion site density on both types of platelets is achieved 
by imposing Nb/Nm = Rb/Rm — 1600; in practice 
Nb = 30400 and Nm — 19, while surface charge densi- 
ties were fixed at cjb ~ 0.09 enm^^ and am — 0.9 enm^^. 
Examples of bentonite pair distribution functions gbbij") 
are shown in Fig. ^ At the two lower clay concentra- 
tions, the position of the first peak scales roughly like 

— 1/3 

, but at the highest density the peak shifts dra- 
matically towards the origin, a behaviour indicative of 
a collapse of nearly parallel pairs of bentonite platelets, 
under the infiuence of a strong effective attraction, in- 
duced by the smaller platelets that appear to " condense" 
on the surfaces of the bentonite discs. Note that all 
secondary minima and maxima have disappeared from 
gpp{r) at the highest density. The corresponding struc- 
ture factors Sbbil) (minus the fixed intramolecular con- 
tribution) are pictured in Fig. ^ together with the results 
obtained from PRISM for the same bentonite densities, 
but ten times lower densities of the smaller platelets. 
The structure factors Sbbiq) at the two lower densities 
show typical liquid-like structure, with a main peak at 

— 1/3 

q = 2tt /Db, where Db ~ Pb ■ Increasing pb or adding 
mixed metal hydroxide shifts the peak toward higher q- 
values, corresponding to a shorter distance between ben- 
tonite platelets. Increasing pm also leads to an enhance- 
ment of the small angle {q 0) value of Sbb{q)] a clear 
central peak appears at intermediate densities, while at 
the highest densities the amplitude of the central peak 
appears to diverge, signalling strong concentration fluc- 
tuations indicative of platelet aggregation, or spinodal in- 
stability reminiscent of the behaviour observed in recent 
computer simulations of mixtures of oppositely charged 
spherical particlesE3. A similar strong small q upturn has 
been observed cxpcrimentalhL for monodisperse clax sus- 
pensions oil the gel phaseEjLj, polyclcctrolyte gelsc£l ami 
mixturesCJ, Poly(styrenesulfonate) ion exchang&j-esinscZl, 
and low ionic strength polyelectrolyte solution^"E3, but 
the exact physical origin of these large-scale fluctuations 
for these systems is still under debate. 
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FIG. 6. Pair correlation function gbb{r) for a mixture of 
large (subscript b) and small, oppositely charged (subscript 
m) platelets for various densities: pb = 0.183 A*m~^, pm = 40 
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FIG. 7. Platelet structure factor for a mixture of large (sub- 
script b) and small, oppositely charged (subscript m) platelets 
for various densities: pb = 0.183 pm~^ , pmi = 4 /im~^, 
Pm2 ~ 40 fj,m~^ (top curves); pb — 0.074 pra~^, pmi ~ 1-6 
pm~^ , pm2 = 16 fj,m~^ (middle curves); pb = 0.037 pra"^, 
pmi — 0.8 //m , pjn2 = 8 /^m (bottom curves). The dashed 
and solid lines represent the calculations for pmi and pm2, re- 
spectively. The upper four curves are shifted. 
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VI. CONCLUSION 

We have illustrated the mechanism of mesoscopic par- 
ticle aggregation induced by oppositely charged polyions 
on two examples of complex two-component dispersions. 
The partial "condensation" of the small particles leads 
to a considerable reduction of the strong Coulomb repul- 
sion between the larger colloidal particles, which may go 
as far as an effective short-range attraction between the 
latter. The attraction, and resulting aggregation give rise 
to a strong enhancement of the long wavelength concen- 
tration fluctuations of the large particles, which is well 
characterized by the small q limit of the corresponding 
partial structure factor. 

While effective attractions have been widely stud- 
ied theoretically and|-b,y simulations in the case of mi- 
croscopic counterionsEJ, PRISM has allowed us to ex- 
tend such studies to cases where the " counterions" are 
themselves complex objects (semiflexible polyelectrolyte 
chains or clay platelets). PRISM turns out to be an 
accurate and flexible tool to tackle highly asymmetric 
complex systems of charged particles of various shapes, 
which are beyond the search of present day simulation 
methodology. The PRISM integral equations provide 
detailed structural information, in the form of partial, 
site-site pair distribution functions and structure factors. 
Thermodynamic properties can then be deteiaiined via 
the energy route, or the compressibility routeQ, but the 
structural information is not sufficient to give direct ac- 
cess to the equation of state via the viral route, and to 
the free energy necessary to determine phase diagrams, 
which requires thermodynamic integration. Work along 
these lines is in progress, as well as further calculations 
based on model B for the coUoid-polyelectrolyte mixture. 
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